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Abstract 

A systematic study of low energy nuclear structure at normal deformation is carried out using 
the Hartree-Fock-Bogoliubov theory extended by the Generator Coordinate Method and mapped 
onto a 5-dimensional collective quadrupole Hamiltonian. Results obtained with the Gogny D1S 
interaction are presented from dripline to dripline for even-even nuclei with proton numbers Z = 10 
to Z = 110 and neutron numbers N < 200. The properties calculated for the ground states are their 
charge radii, 2-particle separation energies, correlation energies, and the intrinsic quadrupole shape 
parameters. For the excited spectroscopy, the observables calculated are the excitation energies 
and quadrupole as well as monopole transition matrix elements. We examine in this work the 
yrast levels up to J = 6, the lowest excited + states, and the two next yrare 2 + states. The 
theory is applicable to more than 90% of the nuclei which have tabulated measurements. We 
assess its accuracy by comparison with experiments on all applicable nuclei where the systematic 
tabulations of the data are available. We find that the predicted radii have an accuracy of 0.6%, 
much better than can be achieved with a smooth phenomenological description. The correlation 
energy obtained from the collective Hamiltonian gives a significant improvement to the accuracy 
of the 2-particle separation energies and to their differences, the 2-particle gaps. Many of the 
properties depend strongly on the intrinsic deformation and we find that the theory is especially 
reliable for strongly deformed nuclei. The distribution of values of the collective structure indicator 
i?42 = E(4f)/E(2f) has a very sharp peak at the value 10/3, in agreement with the existing data. 
On average, the predicted excitation energy and transition strength of the first 2 + excitation are 
12 % and 22% higher than experiment, respectively, with variances of the order of 40-50 %. The 
theory gives a good qualitative account of the range of variation of the excitation energy of the first 
excited + state, but the predicted energies are systematically 50 % high. The calculated yrare 2 + 
states show a clear separation between 7 and (3 excitations, and the energies of the 2 + 7-vibrations 
accord well with experiment. The character of the 0^ state is interpreted as shape coexistence 
or /3-vibrational excitations on the basis of relative quadrupole transition strengths. Bands are 
predicted with the properties of /3-vibrations for many nuclei having R42 values corresponding 
to axial rotors, but the shape coexistence phenomenon is more prevalent. The data set of the 
calculated properties of 1712 even-even nuclei, including spectroscopic properties for 1693 of them, 
are provided in CEA website and EPAPS repository with this article 1]. 
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I. INTRODUCTION 



A present-day goal in nuclear theory is to develop a universal theory of nuclear struc- 
ture, in the sense that it is well-founded in its methodology and is applicable across the 
chart of nuclides. The most promising starting point is self-consistent mean-field, but the 
theory must be extended in some way to treat excitations and nuclear spectroscopy. For 
any candidate theory or methodology, one needs to know its performance on known observ- 
ables to be confident about predictions to unknown nuclei or regions of the nuclear chart. 
It is our purpose to provide and document this information for one particular theory, the 
constrained-Hartree-Fock-Bogoliubov (CHFB) theory together with a mapping to the five- 
dimensional collective Hamiltonian (5DCH), using the Gogny D1S interaction in the nuclear 
Hamiltonian . The results presented here are a major extension of the our study of one 
particular observable in the CHFB+5DCH theory, the low-lying quadrupole excitation J4]. 
Since that work and during the course of the present calculations two new parameteriza- 
tions of the Gogny force have been published 5]. These are aimed at removing systematic 
deviations of calculated binding energies from measurements using the D1S force. However, 
for our purposes, the global assessment of spectroscopy properties, it is more important at 
this stage to benchmark the predictions of a stable and widely tested Hamiltonian. The 
present results will therefore serve as a baseline for comparisons with the next generation 
of similar structure calculations based on the new parameterizations. We mention that the 
present era of global calculations within self-consistent mean field theory using a fixed inter- 

n q 

action started with the calculations of Tajima et al. [6] and Lalazissis et al. [7|. Also, the 
mapping to 5DCH has recently been adapted and applied to the Skyrme and the relativistic 
mean-field Hamiltonians Q , 9] . 

The CHFB+5DCH theory is one of several paths that could be taken to extend mean- 
field theory to describe spectroscopic properties such as excitation energies. Rather than 
constructing a collective Hamiltonian from the CHFB solutions, the constrained wave func- 
tions may be projected on good angular momentum and then used directl y to construct a 
discrete basis Hamiltonian. This technology in now fairly well developed 10|, 111] and has 
recently been applied to the global study of correlation energies 12] and the lowest 2 + exci- 



tations 



13] . In principle, it should give similar results to the CHFB+5DCH theory applying 



the same Hamiltonian. However, in the present implementations the discrete-basis wave 
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functions do not include triaxial deformations and the effective rotational inertias arising 
from the Hamiltonian matrix element are not self-consistent. On the other hand, discrete 
basis methods do not rely on the Gaussian Overlap Approximation which we require for 
our Hamiltonian mapping. There is a third path to introduce dynamics into mean field the- 
ory, the quasi-particle random-phase approximation (QRPA). That has also been applied to 
many nuclei, but so far global surveys have been restricted to spherical nuclei fl^ ]. 

The methodology in the present work has two stages, the CHFB calculations to set up the 
5DCH input, and the solution of the 5DCH equations. For the first part, we perform CHFB 
calculations on a large triaxial grid of quadrupole deformations. These provide a potential 
energy surface and the inertial masses needed for the 5DCH. Positive parity solutions are 
extracted for about 1700 even-even nuclei with proton numbers Z = 10 to Z = 110 and 
neutron numbers N < 200. These calculations span most of the periodic table from drip- 
line to drip-line. One purpose of our work is to establish benchmarks for the accuracy and 
reliability of the theory for energies and properties of the low excitations. We therefore 
make systematic comparisons with experimental data, especially when it is available as 
a tabulation from a published critical review or data repository. The quantities that we 
can easily compare are two-nucleon separation energies and gaps, excitation energies of 
the lowest excited states including yrast spectra up to J = 6, and transition rates of the 
lowest 2 + and excited + states. Since no effective charge is involved, our predictions are 
free of parameters beyond those contained in the Gogny D1S interaction. We hope these 
calculations will be helpful to understand the limitations of the theory and ultimately find 
improved methodologies and Hamiltonians. One important question deals with shell gaps 
and magic numbers. Far from stability, dedicated experiments have shown that the N = 
20, 28 shell gaps experience erosion and that N = 16 may become magic number at the 
oxygen neutron drip-line [ijj]. Other experiments are underway to investigate whether shell 
quenching takes place too in the vicinity of the N = 50 gap, a critical issue for the path 
followed by the r-process of the nucleosynthesis. The predictions for the shell gaps and 
associated observables coming from the Gogny interaction have previously been reported 
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Another purpose of this work is to provide a set of predictions for nuclei to be studied 
in future. The advent of unstable nuclear beam facilities has opened up a new and exciting 
area in exploring the structure of exotic nuclei. Of particular interest are the nuclei near or 
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at the border lines of stability, the proton and neutron drip-lines. These nuclei often display 
properties which are not present in nuclei located in the vicinity of the /9-stability line, and 
questions are raised as to whether the nuclear structure models and effective forces tailored 
over the past seventy years remain valid in the present context. Thus, strong deviations of 
the experimental findings with respect to the benchmarked predicted accuracy would signal 
new phenomena in nuclear structure. 

Several caveats should be mentioned that limit the domain of validity of the theory. A 
basic approximation made here is to require the HFB fields to conserve parity and signature. 
There is no strong evidence that these symmetries are violated in the HFB ground states, but 
there may be some nuclei for which it happens. Another limitation arises from the neglect 
of two- and higher-order- quasiparticle (qp) excitations. The GCM theory (with or without 
the Gaussian overlap approximation) does not include these degrees of freedom which will 
inevitably affect the spectrum at higher excitation energies. Also, the Hamiltonian is an 
adiabatic one, with parameters calculated in the vicinity of zero rotational frequency. This 
affects the reliability of the calculated excitations with high angular momentum. Finally, 
the application of the 5DCH requires further that the overlaps be semilocalized in the (f3, 7) 
quadrupole deformation plane (coherence length be small compared to dimensions of the 
arena in 7) in which the wave function has significant amplitude). Indeed, the mapping 
of the HFB to the collective Hamiltonian is problematic for very rigid nuclei, such as the 
doubly magic ones. For these reasons, we concentrate on the lowest excited states in this 
work in non-doubly magic nuclei, and restricting angular momentum to J < 6. 



II. REMINDER OF FORMALISM AND COMPUTATIONAL IMPLEMENTA- 
TION 

For completeness, we recall here the equations to be solved. The derivation and some 
aspects of the implementation are presented in more detail in Ref. The potential energy 
surface that goes into the 5DCH is obtained from constrained Hartree-Fock-Bogoliubov 
calculations (CHFB) based on the Gogny D1S interaction. The CHFB equations to solve 
are 

5(<l>(q ,q 2 )\H - \ Q - \ 2 Q 2 - \ Z Z - \ N N\<5>(q Q ,q 2 )) = 0. (1) 
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H is the Hamiltonian, and the other terms are linear constraints to obtain particle numbers 
N, Z and quadrupole moments go, 12 according to 

($(9o,g2)|4l$(go,g 2 )) =*, (2) 

($(q ,q 2 )\Z(N)\<f>(q ,q 2 )) = Z(N). 

Here we define the quadrupole operators as Q = 2z 2 — x 2 — y 2 and Q 2 = x 2 — y 2 . The CHFB 
equations are solved for each set of deformations by expanding the single particle states in 
a triaxial harmonic oscillator (HO) basis. 

The triaxial oscillator basis is subject to truncation according to 

(n x + l/2)hw x + (n y + l/2)ftw y + (n g + l/2)fko z < (N + 2)Hlu , 

where (ftcuo) 3 — hjj x hjjj y hjjj z , with hiOiii = x,y,z) as oscillator basis parameters, and rij(z = 
x, y, z) as quantum numbers. The basis is determined as follows. For a nucleus with Z 
protons and N neutrons, the number N is such that A/", the number of single particle states 
in the Hartree-Fock scheme, is eight times the number of levels occupied by the larger among 
the Z or N values. M is a function of N$, fulfilling the empirically established equation 

A/' = 2.1iV 2 + 0.00727V 4 , 

from which N is deduced. In general N Q is not an integer. 

Instead of using the hwi parameters, we adopt in the HFB calculations the parameters 
Hcoq, P and Q, with P= hw x /hw y and Q=hjj x /hjj z . These parameters need be determined to 
define the oscillator basis at each point of the grid. The parameters P and Q are determined 
using formulas based on a liquid drop parametrization of nuclear shape. These formulas 
depend upon (3 and 7 deformations in the constrained HFB calculations, and write as 

P — exp[— x\/3sm7], (3) 
Q = exp(x[^cos7 - ^shi7]), 

where x = (3/ {2(3+ 1). 

hjjQ is obtained through minimization of the HFB energy. This is made for 7= 0° and 
60° at a fixed (3 to take advantage of our axially symmetric HFB code which is running 
much faster than the triaxial code. To get huo values over the triaxial plane, we use the 
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interpolation formula 

WoO0,7) = ~M&7 = 0°) +^ (/3,7 = 60°)] + i[w„(/3,7 = 0°) - w (/3,7 = 60°)] cos(3 7 ). 

For the basis truncation, a fine tuning of P and Q values is performed so as to maximize 
the number of particle states without altering the numbers of oscillator shells in each of the 
three directions as obtained with Eq.Q. Typically, the number of major shells Af ranges 
from 6 to 16 in the present study. 

The Bogoliubov space is restricted by imposing the self-consistent symmetry Ttt2, with 
7r 2 the reflection with respect to the xOz plane, and T the time-reversal symmetry. The 
HFB nuclear states have also been taken invariant under the left-right symmetry 19|. One 
technical point should be mentioned. Since there are many points to calculate, it is important 
to have an efficient algorithm to perform iterative solution of the CHFB equations. From 
the early days, we found it very helpful in this respect to use first order perturbation theory 
to update the linear constraints. During the iterative procedure the obtained mean value 
Qj differs from the imposed value qf 1 , the corrections applied to the Lagrange parameters 



are |20| 



6X i = y £(M i l 1 )- 1 (qf ) -q j ). (4) 

3=0,2 

The moments Ai of the off-diagonal quadrupole operators in the constrained HFB configu- 
rations are defined as 

where /i, v label quasiparticles with destruction operators r\ and energies and E v , respec- 
tively. 

The potential energy surface is then determined from the expectation value of the Hamil- 
tonian, corrected for the one- and two-body center-of-mass energy 

V(q , <fc) = ($(<Zo, <h)\H - ^l<%o, ffe)). (6) 

It is convenient to use the dimensionless deformation parameters (j3, 7) which are defined 
through (3 = \/57r^- fo "t ^ 2 and 7 = arctan a/3— , with r = 1.2 fm. Typically, the 



constrained HFB equations are solved on the domain (0 < (3 < 0.9 ; < 7 < 7r/3) with 
mesh spacings A/3 = 0.05 and A7 = 10°. 
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The final 5DCH is expressed 



3 

H C oll 



\t. J -T-\ E D-^^D^(B mn )-^ + V(a ,a 2 )-AV(a ,a 2 ), (7) 

, Ok ^ ( Ja m <J(l n 

k=l m,n=Daria2 



where we have made another change of deformation parameters from (/3, 7) to ao = /3cos7 
and 0,2 = /3sin7, and where D is the metric 21] . There are 3 rotational inertia and 3 
quadrupole mass parameters in the 5DCH. These are all computed from the local properties 
of the CHFB solutions at the grid points. To calculate the rotational inertia we implement 
additional constraining fields uJk to Eq. (1), where </& is the angular momentum operator 
about the k axis. Calling the new self-consistent solution we calculate the inertias 
as 

Jk = { «' *' g) . (8) 

UJ 

In the limit uj — > 0, this expression is equivalent to the Thouless-Valatin inertia. In practice 
we take u = 0.002 MeV to approximate the limit. The quadrupole mass parameters are 
calculated in the cranking approximation 181 ] . 

BM - - M ' JM , (9) 
tjKq) 2(M t l 1 (q)) 2 

with A4 l jj!(q) the moment defined in Eq.fjSD- 

It is important to mention that the cranking approximation is not self-consistent in the 
sense that the dynamical rearrangement is not taken into account and we should expect 
some deficiencies in the theory result. 

The zero-point energy (ZPE) correction to the potential, Eq.([S]), is associated with the 
nonlocality in the quadrupole coordinates. It is calculated according to the formulas given 
in Refs. (lfibA 

M%(q) 



4 ^ M l lM 

Here the sum runs over the sets = (0, 0), (2, 2), (0, 2) for the vibrational ZPE and 

= (1, 1), (-1, -1), (-2, -2) for the rotational ZPE, following the notation of 0. This 
includes only the part of the ZPE arising from the kinetic energy operator. There is also 
a part due to the potential, which we neglect. This is expected to be small in typical 
situations with shallow minima in the potential energy surface; it might be significant near 
magic numbers where the curvature of the surface is higher. 
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Eigenstates and eigenenergies are obtained as numerical solutions of 

H coll \JM) = E(J)\JM). 



(10) 



The orthonormalized eigenstates \ JM) with angular momentum J and projections M on the 
third axis in the laboratory frame are expanded as 

\JM)=Y J 9 J K^a 2 )\JMK), (11) 

K 

with \ JMK) a superposition of Wigner rotation matrices. The probability P{K) of the 
different K components of the wave function gives a useful indicator of its character. This 
is defined as 

P(K) = J da da 2 \g J K {a^a 2 )\ 2 . (12) 

We refer to Ref. 18| for further numerical details on solving the 5DCH equations. Here 
we use the value m max = 28 for the order parameter in the power expansion of vibrational 
amplitude. This secure a 2% precision on relative energies in collective spectra. We calculate 
radii and quadrupole matrix elements assuming that the coordinate operators are local in 



the collective coordinates 
M- m = r f^m(^) i s calculated as 



211 ] . For example, the matrix element of the quadrupole operator 



(J'M'\M m \JM) = ((2 J + 1)(2J' + 1)) 1/2 x (13) 

J 1 2 j\ I J' 2 J\ f , , 

2 , / daoda2gK,(a ,a2)g K (a ,a2)(§(a ,a 2 )\Mk\$ (a ,a 2 )). 



K,K'k 



—M' m M 



(14) 

This is an approximation, but we have no reason to doubt its reasonableness. 
The correlation energy is defined as 

E CO rr = Eft % p B — E 5 £)CH, (15) 

where E^ l p B is the minimum of the energy at the HFB level, and E 5 dch is the energy of 
the collective ground state obtained from the 5DCH calculations. For nuclei near magic 
numbers, the calculated correlation energy may come out negative, which is unphysical. We 
have kept these nuclei in the accompanying table, but we exclude them when we compare 
the calculated properties with experiment. 
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Also, the accuracy of the calculations will not be as high at the extremes of the nuclear 
chart, due to the incipient shape instability associated with fission, as well as the limita- 
tions of the harmonic oscillator basis for dripline orbitals. In the accompanying tables, we 
include the ground state properties when the calculated deformation is consistent with a 
non vanishing fission barrier in the 7) plane. This will include some nuclei that would 
have vanishing fission barrier when more shape degrees of freedom are permitted. 



III. EXAMPLES 



To show the scope of the theory, we begin with two examples of nuclei that illustrate the 
complexity of nuclear structure that can be addressed with the 5DCH. The first is 76 Kr , 
which is considered as an example of a soft nucleus. The second is 152 Sm, which has a near 
rotational spectrum but is also considered to be a transitional nucleus. 



A. 76 Kr 

We begin with 76 Kr , a nucleus with a complex spectrum of low-lying excitations providing 
evidence for shape coexistence phenomena. A few calculated spectroscopic properties of this 



nucleus were already reported in Ref. 23j. The ground state of 76 Kr is spherical in the HFB 
approximation, but becomes highly deformed in the CHFB+5DCH wave function, with 
mean deformation values of (j3) = 0.33 and (7) = 24°. The variances of the deformations 
are also of interest, namely 



5P = V(P 2 ) ~ (P) 2 ; h = V(l 2 ) - <7> 2 , (16) 

where < 7 2 > and < 7 > are calculated over the sextant < 7 < tt/3. The values of these 
quantities in the 76 Kr ground state are 5f3 = 0.10 and £7 = 13°, suggesting that the nucleus 
is fairly rigid in (3 but with some soft triaxiality. Due to the triaxiality, one does not expect 
to see a rigid rotor spectrum, despite the large deformation. 

In this work, we will examine systematically the Of, 2f, O^ - , 2f, 4f, 2f, 3f, and 6^ excita- 
tions. These are shown for the 76 Kr nucleus together with the additional states that could 



form a 7- vibrational structure in Fig. Q] 23|, |24j . The experimental spectrum is also shown 



in the figure, and one sees that excitation energies are reproduced very well. The calculated 
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FIG. 1: Experimental and theoretical spectra (MeV) and transition strengths (e 2 /m 4 ) of 76 Kr 
showing the excitations that we examine in the present global study. The experimental spectrum, 



on the left, is from Ref. 

band y. Calculated values are those from Ref. 



231 ] as well as from data rep ository for the 3 + and 4 + members of the 7 



23(]. 



excitation energy of the first excited state, the 2f, is only 21% higher than experiment. The 
other states are proportionally even closer: the energies of the yrast 4f and 6f states are 
within 10% of the experimental values. Even the non-yrast excited state energies come out 
well. We shall later examine systematically the 0^~,2j" and 2jj~ excitations; in 76 Kr their 
predicted energies are all within 20% of experiment. For the transition strengths, the pre- 
dicted B(E2; 2f — > Of) is within 20% of the experimental value and the higher transitions 
along the yrast ladder are within 10%. 

The second excited state in the 76 Kr system is the 0^" level at 0.77 MeV. The calculated 
energy is 0.92 MeV, close enough to make a correspondence between the two states. Its mean 
deformation parameters are close to those of the ground state, suggesting a /5-vibrational 
interpretation. The transition rate to the 2f state is large and in very good agreement with 
experiment. The 2% excitation corresponds in excitation energy fairly well to experimentally 
measured state. The 2% wave function has a large probability P(K = 0), suggesting that it 
be placed with the 0^ level as member of the K = excited structure. Its transition strength 
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to the is large and in qualitative accord with experiment. Experiment and calculation for 
the spectroscopic quadrupole moment Q(2g ) are also in accord for both magnitude and sign. 
The sign is opposite to that for Q(2^~), nullifying any interpretation of the K = excited 
structure as /3-vibrational band and giving weight to the interpretation of shape coexistence 
between prolate and oblate band structures. 

From the energetics, the l\ level might be assigned at a two-phonon excitation of the 
ground state. The calculated 2% wave function has a large probability for K = 2 (-P(2) = 
0.77), suggesting that this level instead is the bandhead of a 7-structure. However, the 
experimental data on the transition strengths between the 2^ and the 0^, 2f and 3f states 
are very far from the theoretical predictions. Since transition strengths of 7-vibrations are 
very sensitive to K-b&nd mixing, the disagreement of transition strengths does not rule out 
the 7- vibrational interpretation. For more discussions, see 23]. 



We finally mention the highest excitations, some of which will be beyond the scope of 
our global survey. There is good accord between theory and experiment for the energetics 
of the 6^ state as well as for those for the 3^ and levels which both form a quasi-7 band 
structure on top of the 2^ state. The 6^ and states have strong transitions to the 4+ 
and l\ states, respectively, in both theory and experiment. 

To summarize, the CHFB+5DCH theory provides a very good description of low-lying 
excited states 76 Kr spectrum. While not all aspects are reproduced, many of the energies 
and relative transition strengths are given to good accuracy. The complex spectra of the Kr 
isotopes have often been discussed as a shape coexistence phenomenon, and the theory does 



rather well in describing these features as well as shape transitions in this region [251 ] . 



B. 152 Sm 

The nucleus 152 Sm lies at the start of the deformed lanthanide region of the nuclear chart 
and is considered as landmark in the identification of first-order quantum phase transition 



between spherical and axially deformed nuclei 26.l27l| . Its experimental level scheme is shown 



in the left-hand panel of Fig. [2l taken from Ref. 24| and our calculated level scheme is in the 
right-hand panel. We first note that the yrast band is well reproduced. The l\ excitation 
energy is within 2% of the experimental one, and the experimental ratio of the 4J" to the l\ 
energies is -R42 = 3.0, slightly lower than the rigid axial rotor value 10/3. The 5DCH ratio 
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FIG. 2: Experimental and theoretical spectra (MeV) of Sm. The experimental spectrum is 



based on Ref. 



24] . 



is 3.0, reproducing the slight deviation from rigidity. The yrast spectrum is intermediate 
between that for harmonic vibrators and axial rotors, consistent with predictions from the 
X(5) model designed as a naly tic description of critical point structures in N ~ 90 isotones 
29j. For a review see [30J. 
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We now come to the predictions for the 0^~ excitation and the collective structure built 
upon it. The calculated deformation of that state is (J3) = 0.29, almost the same as the 
ground state deformation, ((3) = 0.30. This suggests an interpretation as a /3-vibration. 
One expects that the fluctuation in j3 would be larger in the vibrational excitation than in 
the ground state; in the harmonic limit, < n\/3 2 \n >= (n + 1/2) /B oq uj, giving a ratio of 
\/3. In fact, the fluctuation in the calculated wave functions is larger for the excited state, 
by a factor 1.58. Thus, the theoretical wave function has the main characteristics to be a 
/3- vibration. 

Comparing experimental and calculated spectra, we first note that the calculated exci- 
tation energy of the 0^ state is quite a bit higher than observed experimentally, by nearly 
40%. As we will see later, this is a common feature of the CHFB+5DCH theory as presently 
implemented. The experimental energy splitting between the l\ and 0^ states is nearly 
identical to that of the ground state band, while the theoretical splitting is larger by 40%. 
This discrepancy is in keeping with that of the X(5) model 30J. We conclude that the 
theory confirms in an approximate manner the existence of a band structure based on the 
0^~ excitation, but in detail deviates from the /^-vibrational limit in the in-band energetics. 

The CHFB+5DCH theory also predicts a 2jj~ excitation and collective structure upon it. 
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TABLE I: Experimental and theoretical E2 transition strengths B(E2; Jj — ► J A (e 2 /m 4 ) of 152 Sm. 



From left to right : expl and exp2 are for B(E2) experimental data from Ref.[31| and 
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341 ] , respectively, and th are for 5DCH calculations. 
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This sequence is interpreted as a quasi 7-vibrational band, with head level energy slightly 
higher than that observed in this nucleus. Our calculated third 2 + state has the same 
average (/?) deformation as the ground state, supporting a vibrational interpretation. If it 
were a true 7-vibration, it should have high probability for the K = 2 component of the 
wave function. This probability is 0.64 compared with 0.002 and 0.35 for the 2+ and 2j 
levels, respectively. Thus, the 2g state has a qualitative character as a 7-vibration but this 
is diluted by other components. This is to be expected for a transitional nucleus such as 
152 Sm. 

Important indicators for structure properties are the strengths for intra- and inter-band 
E2 reduced transition probabilities. B(E2;Ii — > //)'s measured by the Georgia Tech. and 
Yale collaborations are shown in Table 1 together with CHFB+5DCH calculations. As the 
two sets of experimental data display differences, comparison between B(E2) predictions and 
measurements necessarily has a global character. The figure of merit of our theory for 152 Sm 
is as follows : i) the intraband transition strengths have right order of magnitude, especially 
for the ground state band, ii) the transition between 7 and ground state as well as between 
the 7 and (3 bands are too collective, and iii) the f3 to ground state band transitions display 
a mixed character. The theory for the 0% — > 2+ transition strength is about 60% too high. 
Nevertheless we conclude that the 0^~ excitation is a /3-vibration in 152 Sm. 

Our present conclusion is that the structure of the ground state, j3-, and quasi-7 bands is 
globally as the 5DCH theory predicts, but there are probably other components in the wave 



functions, such as 2qp excitations and pairing isomerism 35[ that may have an important 



large effect on the out-of-band transitions. More accurate B(E2) measurements that are 



underway [36|] will be a valuable asset for making definite statements on the predictive 
character of present 5DCH calculations and for disclosing which degrees of freedom might 
be missing in the structure models including the 5DCH one. 

IV. GROUND STATE PROPERTIES 
A. Nuclear shapes 

We begin by displaying in FigJ3]the set of nuclei that we have calculated and included in 
our tables. This comprised all even-even nuclei that are stable with respect to two-particle 
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FIG. 3: (Color online) Chart of nuclides showing ground state deformations. Panel (a): HFB 
minimum; panel (b): expectation value in the 5DCH ground state. The black curve shows the 
beta-stability line. 

emission, and that have positive correlation energies in the CHFB+5DCH theory. Among 
the Z > 96 nuclei, some ones close to the proton drip line have been removed from the chart 
as their inner potential barriers are too low for inhibiting fission decay. For now, we remark 
that the two-particle stability is almost completely determined by the HFB energies. The 
correlation energy E corr contribution (see Eq J15jl ) only changes it for a few nuclei on the 
borders. The next general remark is that the criterion of positive correlation energy affects 
only nuclei at magic numbers. These are visible as the absence of colored circles along some 
of the magic number dotted lines. 

In the figure, the color coding shows the deformation (3 of the ground state, with a ± 
sign according to the value of 7. For the HFB ground states, shown in the left-hand panel, 
one sees the familiar landscape of nuclear shapes, with nuclei near magic numbers having 
small or vanishing deformation (dark and medium green), and two large deformed regions 
located at the lanthanides and actinides. Additional regions of deformation are centered at 
nuclei with (Z, N) = (12, 12), (38, 40), (40, 60), and (60, 80), the heavy nuclei with N ~ 150, 
and the superheavy nuclei with N > 190. It is also apparent that single magic numbers 
do not enforce sphericity. For example, the Sn isotopes are spherical in the region below 
iV ~ 82, become deformed for neutron numbers in the range N ~ 100 — 112, and get back 
to spherical shapes beyond N ~ 114 up to the neutron drip-line. 

The right-hand panel shows the expectation value {(3) for the CHFB+5DCH calculation, 
with the sign determined by the expectation value of 005(3(7)). The 5DCH wave functions 
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have larger deformations on average with fewer nuclei near sphericity. To better see how the 
5DCH changes the deformation properties, we show in Fig. HJhistograms of the distributions 
of (3 and 7. The result for the distribution of (3 in the HFB theory is shown on the upper left- 
hand panel. Among the 1712 nuclei in the calculated data set, roughly 30% are spherical. 
The rest has a broad distribution of deformations peaking at f3 ~ 0.25. Except for the very 
heaviest nuclei, the largest deformation of the survey was found for the nucleus 24 Mg, with 
(3 = 0.54. The lower left-hand panel shows the corresponding distribution of 7. For this plot, 
we restricted the nuclei to those with (3 > 0.1, because 7 is ill-defined in spherical nuclei. 
One sees that the great majority of the nuclei are prolate and axially symmetric, i.e. 7 ~ 0. 
There is also a small peak for oblate shapes, 7 = 60°, comprising about 15% of the deformed 
y of oblate deformations compared to prolate is well-known in mean-field 
. However, it should be mentioned again that our calculated nuclei include 
only those having positive correlation energies. The others are all near magic numbers and 
are likely to be spherical. Turning to the 5DCH (/?) distributions shown in the upper right- 
hand panel of Fig. HI we see essentially all the nuclei become deformed, with deformation 
broadly distributed in the range 0.05 < {(3) < 0.4. The corresponding distribution of axial 
asymmetries (7) in the lower left-hand panel shows that axial symmetry disappears in the 
5DCH wave functions, with average asymmetries going up to 30°. 

Additional information about shape fluctuations is provided by the variances in the de- 
formation parameters, Eq. (Tl6l) . In principle, the value 7 = 30° could arise from a potential 
energy surface that is very soft in the 7 coordinate or from one that has a strong triaxial 
minimum. Fig. [5] shows the distribution of rigidity measures (f3)/5/3 and (7)/^7 for j3 and 7, 
respectively. The /9-rigidity goes to very high values, ((3)/ 5 (3 ~ 10 in the deformed actinides. 
We will find that such high values are present when the nucleus has a well-developed rota- 
tional spectrum. On the other hand, the 7-rigidity is much smaller and is never more than 
~ 3. Without a clear peaking at very large values, it will be problematic to characterize the 
nuclei in terms of the simple models for triaxial shapes. 



nuclei, ihe pauci 
calculations Kj, 37 



B. Radii 



We now examine the predicted charge radii, which we compare with the tabulated ex- 

J. The mean square charge radii r\ are calculated as 



perimental data from Refs. [38 
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FIG. 4: Distributions of f3 and 7 in ground states. Panels a) and c): distribution of the HFB 
minima. Several nuclei at Z ~ 108 have minima at (5 ~ 0.75 and are not shown. Panels b) and d): 
distributions of and (7) in the CHFB+5DCH ground state. The distribution of 7 in the lower 
left-hand figure includes only nuclei with nonspherical minima. The histogram on lower right-hand 
includes nuclei having spherical HFB minima as well. Units for 7 and (7) are degrees. 



1 f o t3 / \ 2 ^ 2 

- / r d rn p (r) +r p + —r n 



(17) 



where n p (r) is the point-proton density, and r p = 0.63 fm 2 and r\ = —0.12 fm 2 are the rms 



proton and neutron charge radii, respective 
r 2 m = 3h 2 /2muA fm 2 (see Eq.(4.3) in Ref. 



y. The center-of-mass correction is computed as 



40j), with u = 1.85 + 35.5/A 1 / 3 MeV. We show 



in Fig. [6] the comparison of calculated and experimental charge radii, plotted as the relative 
error 

e = r*Vrf-l. (18) 
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FIG. 5: Distribution of rigidity parameters, namely {(3) /5(3 in panel a) and (7)/#7 in panel b). 

The upper and lower panels show the HFB and the CHFB+5DCH results, respectively, with 
lines connecting nuclei in isotopic chains. We see that the theory is remarkably accurate at 
the HFB level, and the CHFB+5DCH hardly changes the predictions. Among the heaviest 
nuclei, we find that the U isotopes are reproduced very well. The theory seems to be high 
for the Cm isotopes, but it should be noted that these radii were based on systematics in 



the absence of any direct measurement 



4l|. 



Nucleus-to-nucleus variations in radii can be attributed to deformation changes 



421 as 



well as other nuclear structure effects 



43 



44 



45| . The effects of deformation can be easily 



seen in individual isotopic chains. An example is the Sr isotopic chain, shown in Fig. [71 
Experimentally, one sees a slight decrease in the radius from N = 40 to the N = 50 magic 
number, followed by a much steeper increase in radii as more neutrons are added. The HFB 
minima are spherical below N = 50 and deformations increase to very large values at the 
heaviest isotopes in the figure. That results in almost monotone increase in radius from the 
lightest to the heaviest isotopes. 

Turning to the CHFB+5DCH results, we find that the main effect is in the lighter nuclei, 
and it is to increase the charge radius. This is to be expected, since deformations increase 
the radius and the average deformations are systematically larger in the CHFB+5DCH. 
The largest increase, by 4%, is in the nucleus 30 Si. Here the HFB minimum is spherical, 
while the CHFB+5DCH ground state has a mean deformation ((3) = 0.48. Returning to 
the Sr isotopic chain, the correlations associated with the CHFB+5DCH bring the theory 
in very good overall agreement with data. This comes about from two effects. In the very 
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FIG. 6: Charge radii. Plotted are the relative errors, Eq. (|18p . with isotopic chains connected by 
lines. Panels a) and b) show the 
Experimental data is from Refs. 



lines. Panels a) and b) show the results of the HFB and the CHFB+5DCH theories, respectively. 



M, |4fi|, |47( (see 



39l |) and includes 313 nuclei. 



light isotopes, the CHFB+5DCH predicts large deformations instead of the spherical shape 
of the HFB minimum, increasing the radii. On the other end of the isotopic chain the 
nuclei are also deformed, but the average deformation in the CHFB+5DCH wave functions 
~ 0.3 - 0.35) is less than in the HFB minima (J3 ~ 0.45). 
Table [II] shows the performance of the theory, using as a quantitative measure the rms 
dispersion a about the mean e, a = ((e - e) 2 ) 1 / 2 . Both HFB and the CHFB+5DCH treat- 
ments, suitably renormalized, are accurate to 0.6%. For a comparison, the 2-parameter 

taking r c = TqA 1 ^ + r^ -1 / 3 fm is shown in the third row of the 



'Finite surface" model 



Table. Here the error is about twice as large. 



C. Correlation energies 

A key observable that theory should describe is nuclear masses or equivalently their 
binding energies. We shall consider the binding energy to be composed of two terms, the 
binding energy of the mean-field minimum calculated in an unconstrained (with respect to 
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FIG. 7: Charge radii r c of Sr isotopes. Experimental: circles joined by solid line; HFB: triangles 
joined by dot-dashed line; CHFB+5DCH: squares joined by dotted line. 



TABLE II: Comparison of calculated charge radii with experiment: e is the mean of e (see Eq. 
(|18p ): a is its rms dispersion about the average. 313 nuclear radii were included in the comparison 
as in Fig. [6l In the column "HFB (new)" we use the modern value r p = 0.875 fm for the proton 



charge radius 



Theory 


e a 


HFB 
HFB (new) 
CHFB+5DCH 
Finite Surface 


0.001 0.006 
0.005 0.007 
0.006 0.007 
0.0000 0.012 



shape) HFB calculation, and the correlation energy associated with the spread of the wave 
function over the quadrupole shape degrees of freedom. For an orientation, we show in Fig . 
[8] these two contributions and their sum, displayed as difference between experimental 49] 
and theoretical energies (i.e. residuals). One can see that the shell effects at N = 82 and 
126 are too large in the HFB theory, and the correlation energies vary in a way to reduce 
the shell effects to a level closer to that needed. 

The overall performance of the theory with respect to masses depends extremely sensi- 
tively on the parameters of the functional, and any useful theoretical mass table requires 
that the force parameters be refitted. This has been recently carried out for the Gogny 
DIN and DIM parametrizations [^j. However, in the present study we will keep the original 
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FIG. 8: Panel a): residuals of the HFB binding energies with respect to experiment, plotted as a 
function of neutron number with nuclei of the same Z connected by lines. Panel b): Correlation 
energy, Eq. (|15p . Panel c): residuals including correlation energies. Experimental data is from 



Ref. 
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D1S interaction and evaluate the performance with respect to differential quantities, which 
are much less sensitive to the precise parameters of the interaction. The first quantity we 
examine is the two-nucleon separation energy defined as 

S 2n {N, Z) = E(N -2,Z)- E(N, Z), 2n separation energy, 
S 2p {N, Z) = E(N, Z - 2) - E(N, Z), 2p separation energy. 
In the left- and right-hand panels of FigJS] we show the two-nucleon separation energies S% n 
and S 2p for the HFB and the CHFB+5DCH calculations, presenting the calculations in a 



similar way as was done in Ref. 50J. The shell gaps are quite obvious, and one can see that 
they are reduced in the CHFB+5DCH theory. The available experimental data is shown 
on the bottom panels. One can see that the shell gap varies with the number of nucleons 
of opposite isospin. In particular, it is observed in the right-hand panel for the proton 
separation energies that the Z = 50 and Z = 82 shell gaps disappear at high neutron excess. 
As mentioned earlier, the ground states become deformed in these neutron-rich nuclei. 

In the left-hand panel for S 2n one also sees a gradual opening of the N = 162 spherical 
shell gap for proton numbers Z > 96. This gap 2.5 MeV wide for Z = 110 should increase 
stability of superheavy elements (SHEs). Our predictions are consistent with those based on 



calculated shell correction energies 51], and with the observation of a minimum in alpha 



decay energies of SHEs at N — 162 (for a review see 52|). Impact of this neutron gap on 
calculated S 2p values is also seen for iV ~ 162 in the right-hand panel. Finally we note that 
Interacting Boson Model calculations are also supporting evidence for a neutron spherical 
gap in close vicinity of N = 162 531 ] . 



Another phenomenon is the enhancement of the gap near doubly magic nuclei. This 



phenomenon, called "mutually enhanced magicity" 54j, is reproduced much better by the 
CHFB+5DCH theory than by the HFB. The best example is the Z = 82 gap of the S 2p 
systematics in the bottom right-hand panel, which becomes larger near N = 126. Unfor- 
tunately, the Gaussian Overlap Approximation does not permit us to calculate the doubly 
magic nuclei. 

In Table II I II we show the rms residuals of the calculated separation energies with respect 
to experiment. The experimental data is from Ref. j^, including only nuclei whose binding 
energies are given with experimental error of less than 200 keV. As already mentioned, our 
theory only includes nuclei whose correlation energy is positive. This excludes only about 
10% of the nuclei in the experimental data set. The number of nuclei in the comparison 
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FIG. 9: (Color online) Panel (a) : two-neutron separation energies; panel (b): two-proton separa- 
tion energies. The three rows show the HFB theory, the CHFB+5DCH theory, and experiment, 



respectively. Experimental data is from Ref. 



49] 



is given on the first line of the table. The first comparison, with the HFB energies, shows 
rms residuals of slightly less than 1 MeV for both separation energies and gaps. The per- 
formance here is slightly better than was found in the survey based on the Skyrme energy 
functional Sly4, reported in Ref. Jjj]. The bottom line of the table shows the energies of the 
full CHFB+5DCH theory, i.e. with the correlation energy included. The improvement is 
about 25%. This is surprisingly comparable to the results found in Ref. [ijj], despite that 
correlation energy was calculated in a completely different way. 

We also carried out the statistics on the two-nucleon gaps. This quantity is denned by 
the next higher order difference, 

5 2n (N, Z) = S 2n (N + 2, Z)- S 2n (N, Z), 2n gap, 
5 2p (N, Z) = S 2p (N, Z + 2)- S 2p (N, Z), 2p gap. 
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TABLE III: 2-nucleon separation energies and gaps. Sizes of the compared data sets are given on 
the first line. Rms residuals with respect to experiment are given on the third and fourth lines, for 
the HFB and CHFB+5DCH theories, respectively. Energies are in MeV. 





5*2n S2p &2n hp 


Size theory 
exp. 


455 433 396 358 
492 467 444 392 


Theory HFB 

CHFB+5DCH 


1.00 0.91 1.06 0.98 
0.72 0.71 0.68 0.61 



As a particular example, there has been much discussion of evolution of the Z = 28 gap for 
high neutron numbers. We find that the CHFB+5DCH energies are below the HFB values, 
thus weakening any shell effect at Z = 28. There is a peaking at A^ = 28 that could be 
attributed to "mutually enhanced magicity" or to an Z = N symmetry effect, the "Wigner 
energy". Experimentally, there is a slight peaking in the gap at A^ = 40, but we find that it 
is smooth in the CHFB+5DCH theory. 

The overall statistics for the performance of the theories with respect to two-nucleon gaps 
are also shown in Table IHIi The results are somewhat better than those for the separation 
energies. 



V. YRAST SPECTRUM 



In this section we report the predictions for the lowest excitations of angular momentum 
J = 2,4 and 6. For the quantitative measure of the global performance of the theory, we 
will use the same figures of merit as in Ref. 4J. Because the quantities span a large range 
values, we examine the statistics of the logarithmic ratio of theory to experiment, namely 

R x = \og(xth/x exp ), (21) 

for a quantity x. We present its average over the data set R x as well as the dispersion about 
the average, 

a x = ((R x -R x ) 2 } 1/2 . (22) 

The results for the properties we can compare with tabulated experimental data are discussed 
individually below and in Sec. PVTl and are summarized in Table HVl in Sec. IVHI 
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A. The first 2 + excitation 



The first physical property we examine is the fraction of the energy-weighted sum rule 
(EWSR) contained in the 2f excitation. That quantity is governed more by the inertial and 
mass properties of the CHFB+5DCH than by the topology of the potential energy surface. 



The sum rule fraction is often expressed with respect to Lane's isoscalar sum rule [55 ] 



where m is the nucleon mass and (r ) is the mean 



25 h 2 

S(I) = E(2i)B(E2; 0+ - 2+) = ^(-)^ 2 >, (23) 

i 

square mass radius. It is also common 
to make the approximation (r 2 ) = 1.2 2 A 2 ^ 3 fm 2 56] but we shall rather use our calculated 
mass radius. The charged part of the isoscalar sum rule is derived from S(I) assuming that 
the charge current and mass current are proportional 

S(II) = S(I)(jy. (24) 

The fraction s(X) of the sum rules carried by the 2f excitations is calculated as 

s(X) = E(2f)B(E2; 0+ -> 2f)/e 2 S(X), with X = I, II. (25) 

Histograms of s(I) and s(II) are shown on the left-hand panel of Fig. [101 Individually, 
the excitation energies and transition strengths vary over several orders of magnitude. But 
their product scaled by s(X) compresses the rms variation down to about a factor of 2. 
This may be seen in the histograms of s(J) and s(II) shown in the left-hand panel of Fig. 
[TU1 The fraction of strength in each sum rule is about 1.5% for S(I) and 10% for S(II). 
As well known, most of the strength is carried by the giant quadrupole resonance. One 
can also see from the histograms that the scaling with S(II) produces more compressed 
distribution than scaling with S(I). A scatter plot of the theory versus the experimental 
values of S = E(2f)B(E2; Of —* 2f) is shown on the right-hand panel of Fig. [TUJ There 
is a concentration of points on the diagonal that show very good agreement; these mostly 
correspond to strongly deformed nuclei. Overall, the theory somewhat overestimates the 
fraction of the EWSR carried by the 2f state. 

We now turn to the comparison of excitation energies and transition strengths with 
experiment. The results were reported already in Ref. 4|, but we since discovered that the 
code we had been using to solve the 5DCH did not have the desired precision for smallest 
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FIG. 10: Panel a): distribution of sum rule fraction s(X), Eq. ([25]) . in the CHFB+5DCH theory, 
for the 1609 calculated nuclei. Solid and dashed lines show the fraction of the S(I) and S(II) sum 
rules, respectively. Panel b): calculated S = E(2f)B(E2\ 0j~ — > 2f) versus experimental, for 311 



nuclei. Experimental data are from Refs. 



56. 



571 ] . S values are in MeV-e 2 6 2 units. 



excitation energies. In the present work we report recalculated energies using a more accurate 



code described in Refs. 



18 



The comparison of experiment 56|, |57[ and the calculation 
is shown in the left-hand panel of Fig. [TTJ The points at the lower left correspond to 
the deformed lanthanides and actinides, and one sees that the theory does very well there. 
Right-hand panel of Fig. [11] shows a similar comparison for the B(E2; Of -> 2f) transition 
strength. The points on the upper right side of the figure correspond to the very deformed 
actinide nuclei. Again, the theory is seen to be remarkably accurate under the conditions of 
a large static deformation. The global performance figures of merit for the 2+ energy and 
the B(E2; Of — > 2±) strength are given on the first two lines of Table HVl posted in Sec. IVIII 



B. The first 4 + excitation and R42 

An important signature of the character of the excitation spectrum is the relationship of 
the lowest 4 + excitation and the 2f below it. A very useful indicator is the ratio of the two 
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FIG. 11: Panel a): theoretical 2^ excitation energies of 537 even-even nuclei as a function of their 
experimental values. Panel b): theoretical B(E2;0f — > 2^) transition strengths of 320 even-even 
nuclei as a function of their experimental values. Several cases showing deviations are labeled by 
the nucleus. Experimental data is from Refs. jo^l. IotI]. 



excitation energies E(J£), 



R 
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(26) 



The R42 indicator has been much used, particularly in discussing complex spectra. The 
value R42 = 10/3 is characteristic of an axial rotor, R42 = 2 of a vibrator, and R42 = 5/2 
of a 7- unstable rotor or the 0(6) algebraic model 59]. In Fig. [E2] we display histograms 
of the experimental and theoretical ratios side by side. One sees a very narrow peak at 
10/3, showing that one can make a nearly unambiguous assignment of axial rotors. In 
the algebraic models the three simple limits mentioned above represent extremes in the 
parameter space of the models, and it is interesting to see which ones are favored in the 
global systematics. While the axial rotor is clearly special, neither the harmonic vibrator 
nor the 7-unstable rotor shows a corresponding accumulation in the experimental data. The 
CHFB+5DCH theory, on the other hand, does show a second peak just below the 7-unstable 
value, -R42 = 5/2. 

It is interesting to see how well the physical structure indicator i? 42 correlates with the 
intrinsic shape properties of the CHFB+5DCH wave functions. Let us first examine the 
relationship between i? 42 and mean deformation (/?). This is shown in the left-hand panel 
of FigJT3j One can see that the value of ((3) by itself does not determine whether the 
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FIG. 12: Panel a): histogram of experimental ratios, Eq. 



with data from Ref. 



, for 501 even-even nuclei, 



2jyj. Panel b): histogram of calculated R42 ratios for 1609 even-even nuclei 



calculated in the CHFB+5DCH theory. 

yrast spectrum has a rotational character. The R42 has the rotational value for in the 
range 0.2 — 0.45, but nuclei with nonrotational spectra are common with (j3) values up to 
0.3. In fact the highest value of (j3) in our calculations is found for a nucleus ( 26 Mg) for 
which R42 = 2.4, both theoretically and experimentally. Evidently, what is needed as well 
to determine the rotational properties is a measure of the rigidity of the shape. For that 
purpose, we use the /5-softness parameter, rp = Sj3/{/3). The i? 42 values are plotted with 
respect to rp in the right-hand panel of Fig. [13j As may be seen from the figure, this provides 
a much better separation between the rotational and nonrotational spectra. Effectively, the 
/5-softness parameter should be less than 0.2 for a rotational spectrum. 

We turn to the performance of the theory of the 4f level, comparing energies to exper- 
imental data. Of the 484 nuclei with tabulated experimental energies 24], 480 meet the 
criteria to be included in our theoretical data base. Left-hand panel of Fig. [TH shows the 
comparison of the theory to experiment as a scatter plot for i? 42 . For most nuclei the R42 
values in both measurements and calculations fall between i? 42 = 2 and R42 = 10/3 limits 
of the vibrational and rotational models, respectively. Values of R42 less than one are cer- 
tainly possible when the spectrum is dominated by two-quasiparticle excitations, which is 
common near magic numbers. Statistical performance for the data set is given in Table IIVL 
The average value of the ratio R42 comes out very well, only 3% higher than the measured 
average. The dispersion about the mean is also quite good, better than the predictions for 
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FIG. 13: (Color online) The ratio R42 as a function of the mean deformation {(3) (panel (a)) and 
the /3-softness parameter 5(3/ (0) (panel (b)) for the calculated nuclei in their ground states. Dotted 
lines show the rigid rotor value -R42 = 10/3. 
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FIG. 14: Panel a): R42 comparison of theory and experiment for 480 nuclei. Arrows indicate the 



rigid rotor value -R42 = 10/3. Experimental data is from Ref. 
for CHFB+5DCH wave functions of 1693 nuclei. 



2J. Panel b): distribution of R± 2 {B) 
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the absolute energies of the 2f excitations. 

In the subsections below, we will analyze the properties of other excitations with respect 
to the rotational character of the ground state. Since the R42 measure is very clear, we shall 
make much use of it to examine the connections. 

It is of interest to examine the transition strengths B(E2; 4+ — > 2+), even in the absence 
of a critical review and evaluated tabulation of the experimental data. To interpret this 
quantity we take the ratio to the 2f — > 0~t transition, defining R^B) as 

B{E2-Aj^2j) 
Ra2[B) ~ B(E2-2i^0lY 

Two anchor points to interpret R^B) are the axial rotor model for which R^iB) = 10/7, 
and the harmonic vibrator model for which R^B) = 2. The distribution of calculated 
values is shown as a histogram in the right-hand panel of Fig. [TH There is a peak at the 
axial rotor value, but no peak at the vibrator value or anywhere else. The calculated i? 42 (i?) 
range from 1.43 for the nucleus 240 Cm to 5.7 for the nucleus 180 Pb. There are no calculated 
nuclei with R i2 {B) smaller than the axial rotor value. 

C. The first 6 + excitation 

The last excitation we shall examine in the yrast spectrum is the 6|" level. If the 0^, 2^, 
and 4f levels form a band with energies close to the axial rotor limit, the 6f state is also 
part of the band in the vast majority of cases. Deviations of its energy from the rotational 
limit can also be extrapolated from the R42 values using the Mailman systematics 60, |6l| , 
namely the empirical correlation of the ratios Rq2 = E(6f)/E(2f) and R^- The correlation 
associated with the CHFB+5DCH energies is shown in Fig. [15| left-hand panel, based on 
theoretical energies from 1609 nuclei. The scatter plot follows a line from about (R42, Rw) = 
(1.5,2.0) to the value (10/3,7) corresponding to the rotational limit. The plot shows an 
accumulation of points at the axial rotor limit, as well as a somewhat broader peaking near 
(2.3,3.7). For orientation, the positions of the 7-unstable limit and the harmonic vibrator 
limit are for (5/2,4.5) and (2,2), respectively. The experimental scatter plot of Rq2 vs. 
-R42 is shown in the right-hand panel of Fig. [151 It shows data for 458 nuclei, obtained 



from the Brookhaven database 



24J. We also show in the middle panel the calculated nuclei 



corresponding to the experimentally known ones. The experimental points form a line very 
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FIG. 15: Ratios of yrast excitation energies, Rq2i &s & 

function of Panel (a): CHFB+5DC 



theory for 1609 nuclei. Panel (c): experimental ratios for 456 nuclei, with data from Ref. 
Panel (b): theoretical values for the nuclei shown in the right-hand panel. 



much like the one seen in the theory plot. The correlation is also very narrow for the upper 
half of the line, but it becomes broader at lower values of -R42 and R$2- The experimental plot 
extends to lower values than we find in the theory. One possible explanation is the neglect 
of two-quasiparticle configurations in the theoretical wave functions. Such configurations 
can produce high angular momentum at relatively little energy cost, and therefore can give 
values of and Rq2 close to 1. Also, the nuclei with such low i? 62 values may have failed 
our criteria to keep in the theoretical database. The global figures of merit of the observable 
Rq2 are reported in Table [TV] The reliability of the theory is quite high, although it does 
not do as well as for the lower 2 + and 4 + yrast excitations. 
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VI. NON-YRAST EXCITATIONS 



This section will examine in some detail the physical properties of O2 and l\ excited states 
but is also concerned with their possible role as head levels of collective bands, traditionally 
referred to as (3- and 7-vibrational bands, respectively. In order to do this we will also need 
to consider the 2jj~ excitation, which can very often be considered as part of a K = excited 
band. The specific indicators we will examine in this context are the excitation energy 
with respect to band head, the in-band transition rate compared to that for the ground 
state band, and the relative out-of-band transition matrix elements. Unfortunately, the 
data tabulations do not exist to make a systematic comparison to experiment. However, the 
band character has been much discussed in the rare-earth region, and we can compare some 
out-of-band rates there. As the 2jj~ levels are almost systematically members of the excited 
K = bands, an alternative to /^-vibrational band interpretation is suggested, namely that 
of coexisting band structure inside nuclei. A detailed discussion of the 7-degree of freedom 
and associated collective band excitations is deferred to a later publication. 

A. The 2^ excitation 

The lowest non-yrast excitation typically has J = 2, and it is often interpreted in the 
collective model as a shape excitation in the 7 degree of freedom. A theoretical indicator 
for that character is the A"-content of the wave function. This is shown visually in Fig. [IB] 
indicating the probability P(K) by the coloring of the nuclides. Apart from nuclei close to 
magic numbers, the vast majority of second 2 + states have P(K = 2) > 0.75 and can be 
considered as 7- vibrations. In the upper right corner of Fig. [16] is a domain without coloring. 
In this narrow mass region the inner potential barrier is not high enough to sustain excited 
states, and the nuclei go to fission. For a more quantitative view of the K distributions we 
show them by histograms in Fig. [T7] for the second and third J = 2 states in the spectrum. 
For the l\ state (left-hand panel), there is a sharp peak close to P(K = 2) = 1, together 
with a broader distribution of lower probabilities. For the most part, the nuclei within the 
sharp peak have R42 close to the axial rotor value. Thus, for these nuclei we have a clear 
identification of the 2^ level as a 7-excitation. 

The plot for the 2g state in the right hand panel, shows that this level may be viewed 
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FIG. 16: (Color online) Chart of the computed nuclei showing the probability of K = 2 component 
in the wave function of the 2^ levels. The 2^" states with more than 75 % of K = 2 components in 
wave functions are considered as 7-vibrations. The black curve shows the beta-stability line. 
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FIG. 17: Probability of K = 2 component in the wave functions of the second (panel a) and third 
(panel b) excited 2 + states. 

as a /3-excitation in many nuclei. Here the strong peak is at P(K = 2) = 0. Interestingly, 
there are a few nuclei for which the roles of the second and third state are reversed, as can 
be seen in Fig [161 It happens that our example 152 Sm in Sec. Ill is of this kind. In the 
discussion below, we will designate the second or third 2 + state with the larger P(K = 2) 
the 2+ level, if P(K = 2) > 0.75 for all nuclei with i? 42 > 2.3 even though we know that 
7— vibration is a designation specific to well deformed nuclei. 

The systematics of the 5DCH 2+ energies are shown in Fig. [18] (open circles) as a function 
of neutron number. The distribution of excitation energies displays sharp structures with 
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FIG. 18: (Color online) Excitation energy of the CHFB+5DCH 2+ level with P(K = 2) > 0.75 as 
a function of neutron number N. Open symbols in blue and red colors indicate the 2^" and the 2g~ 
levels for nuclei with R42 > 2.3, respectively. These levels are defined as 2+ excitations. Dots are 
for nuclei with < 2.3. 

maxima near N> 50 magic numbers, for which i? 42 < 2.3. Minima are found, as expected, 
half way between major closed shells and they reach very low values (E x ~ 200 keV) in 
heavy nuclei with Z> 98. The performance of the CHFB+5DCH on the energies of the 
2^" levels is shown in Fig. [19] through comparing the calculations to the evaluated data for 



354 nuclei 



24j. The theory clearly reproduces the variation of the experimental energies, 



which range over more than an order of magnitude. The colored symbols in this figure 
are for 1\ states identified as 2+ levels in the CHFB+5DCH calculations. On average the 
theoretical energies are somewhat high. The figures of merit, given in Table HVl show that 
the theoretical energies average about 25% higher than the experimental ones. Interestingly, 
the variance a is smaller than that for the 2^ excitations. 
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FIG. 19: (Color online) Excitation energy of the second J = 2 excitation, comparing 352 nuclei. 



Experimental data is from Ref. 



24 ] . The 2+ levels are marked with red color. 



B. The Of excitation 



In the framework of the CHFB+5DCH theory, the 0^" excitation can arise in several ways: 
as a /3-vibration, as a coexisting state of very different shape, or something in between. To 
be a /3-vibration, the excitation should have nearly the same (/3) as the ground state, but a 
larger dispersion 5(3. The excitation in this limit is very dependent on the calculated mass 
parameters for the j3 degree of freedom, which have been calculated using the Inglis-Belyaev 
formula. This treatment has known deficiencies and we expect that the predicted excitation 
energies would be somewhat lower if the Thouless-Valatin prescription were used. The other 
likely structure for the 0^ excitation arises from the coexistence of vastly different deforma- 
tion at nearly the same energy. The latter mechanism is prominent in light doubly-magic 
nuclei 16 2| . |63| and also in the actinides where the superdeformations occur at low excitation 
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66[ . A phenomenological signature of coexistence would be a low excitation energy. 
In fact there are a number of known nuclei for which the level is the first excited state, 
but we do not find such a low excitation energy in our calculations. However this observation 
does not at all mean that the present theory is not able to provide reliable predictions for 
nuclei where shape coexistence and shape transition are present and characterized by many 
measurements. Such features are well described by our theory for the neutron-deficient Kr 
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FIG. 20: Excitation energy of the 0^ state compared with experiment 24[] . 

isotopes 
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67j |. That the present theory does not predict 0^ state as first excited state 
in a nucleus obviously means that degrees of freedom other than collective quadrupole ones 
are at play and cannot be ignored. 

The theoretical and experimental excitation energies of the excitation are compared 
in Fi g. [201 The experimental data set of 332 nuclei was obtained from the Brookhaven data 
base [24|. Of the 332 tabulated nuclei, 317 are in the CHFB+5DCH calculated nuclei and 
are shown in the Figure. One sees that the theory reproduces the overall variation over one 
order of magnitude, but that the calculations are systematically too high. The figures of 
merit for the performance of the theory are given in Table HVl The average Re is given by 
Re = 0.38 and corresponds to predicted energies that are too high by ~ 50%. The rms 
fluctuation about renormalized theoretical energies is given by oe in the table. Its value, 
0.30, corresponds to a fluctuation +35% — 25% in the error. 



We now examine the 0^~ energy as a function of deformation, following the work of Chou 



et al. [68|. These authors observed that there is a strong empirical correlation between 
the ratio i?02 = ^(0^)/ E(2^) and R42 measurements as shown in the right-hand panel of 
Fig. [2TJ The left-hand panel shows a scatter plot of these quantities for the CHFB+5DCH 
calculations. Indeed, one sees a strong correlation between the two ratios. The R 02 ratio is 
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FIG. 21: (Color online) Panel (a): ratio E{0Z)/E(2f) as a function of -R42 for the calculated nuclei. 



Panel (c): experimental data from Ref. 



24(. Panel (b): CHFB+5DCH values for the nuclei shown 



in the right-hand panel. Color code is for mean ground state deformation. 



flat with a value in the range 1-5 until R42 approaches the axial rotor value, and then it can 
become very large. However, the width of the curve at fixed R42 is too broad to use this plot 
in a predictive way. The overall shape of the curve comes mostly from the variation of the 
denominator in the ratio _E'(0 2 f )/-E'(2+) for both measurements and calculations. We show 
in Fig. [22] a plot of -E'(O^) itself versus R42, which is also informative. Here one sees three 
categories of nuclei. At one extreme are the spherical nuclei, having R42 ~ 2, which have the 
highest excitations for the 0^ levels, typically in the range of 2-4 MeV. At the other extreme 
are the axial rotor nuclei with energies in the range of 1-2 MeV. Interestingly, the nuclei 
in between favor lower energies, in the range of 0.5-1.5 MeV. Those nuclei are likely to be 
soft ones, and that would be reflected in both the excitation energy of the O2" levels and the 
range of values of _R 42 . 
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FIG. 22: (Color online) CHFB+5DCH excitation energy (MeV) of the 0^ state as a function of the 
ratio i?42- Color code is for mean ground state deformation. The vertical line (red color) indicates 
the rotational limit R42 = 10/3. 

1. Criteria for the occurrence of (3 -vibration 

As said above, the Of excitation in nuclei can arise as a /3-vibration, as a coexisting 
level with deformation different from that for lower-lying states, or something intermediate 
between these two extreme structures. Here we would like to place the above discussion in 
a broader perspective. 

First we consider relationships between quadrupole transition matrix elements calculated 
in our theory for the 2^" — ^ 0^~, 2j^" — ^ 2^", and 2f — > Of transitions. If the spectrum truly 
exhibits a /^-vibrational band, the quadrupole transitions between it and the ground state 
should be governed by a single parameter, the matrix element of the quadrupole operator 
between the two intrinsic states. Under these circumstances the spectroscopic transition 



matrix elements are related by Clebsch-Gordan coefficients, cf. 
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, Eq. 4-219]: 



(pj p \\M(E2)\\gJ g ) = {2J g + lfl\j g 020\J p 0){(3\M{E2)\g). 



(27) 



Here the ground-state and the /^-vibration bands are labeled by g and /?, respectively. We 
examine now the three CHFB+5DCH cross-band transitions, (J/3, J g ) = (0,2), (2,0), (2,2), 
to see how well Eq. (1271) is satisfied. According to the model, the magnitudes |Mj /3j j g | should 
satisfy 



I M, 



02 1 



|M 



20 1 




M 



22 



(28) 



To display the deviations of the computed matrix elements from these conditions, we take 
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FIG. 23: Crossover matrix elements. Relative magnitudes of the three quantities 
IM02I, IM20I) \J~^q\M22\ are shown by distances to the sides of the triangle. The vertexes of the 
triangle correspond to the case where only the labeled transition is nonzero. 

the ratio of the three quantities IM02I, | A/20 1 5 y^|M22| to their total. The fractions are 
plotted in Fig. [23] as points within a triangle, the fraction given by the distance to a side 
of the triangle. We see that there is a concentration of points at the center point of the 
triangle; of the 1707 calculated nuclei, 398 have values of the relative matrix elements within 
15% of equality. The distribution of these nuclei in Z and N is shown in FigJ2H One sees 
four regions where the condition is well satisfied, including the strongly deformed rare earths 
and actinides. We conclude that the CHFB+5DCH theory predicts that /3-vibrational bands 
should be quite common, taking as a criterion that eq.(26) be approximately satisfied. There 
is also a concentration of points at the upper apex of the triangle in Fig[23j For these nuclei, 
the 0^ — > 2f matrix element is much larger than the two matrix elements involving the 2g~ 
excitation, suggesting that the excitations behave more like independent phonons. 

2. — ► 2f transition 

We now focus on the experimental situation with respect to the 0^ — > 2± transition 
strength. This is an important observable in an ongoing controversy about the existence of 
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FIG. 24: Chart of nuclei (full circles) in the vicinity of the center point of the triangle shown in 
Fig. [23l The continuous curve is for the ^-stability line, and dots are for nuclei between driplines 
as shown in the left-hand panel of Fig. [3l 



/3-bands in deformed nuclei 



m 



In TOf] , it is concluded that the observed (3- to ground state -band transitions are orders of 
magnitude weaker than those predicted by collective models for deformed rare-earth nuclei, 
except for very few. To assess the performance of the CHFB+5DCH at least in this limited 
region, we have compared the CHFB+5DCH calculations of the B(E2; — > 2f) strengths 



with the experimental data on the 9 nuclei compiled in Ref. [70J. For all of these nuclei, 
the CHFB+5DCH theory predicts a /3-vibrational band that satisfies the criteria discussed 
in the previous subsection. For 4 of the nuclei ( 152 ' 154 Sm, 154 Gd, and 168 Yb), the calculated 
B(E2)'s are of the same order as the experimental ones, but somewhat higher by up to a 
factor of 2 or so. However, for the remaining 5 nuclei ( 158 Gd, 166 > 168 Er, and 172 ' 174 Yb), the 
experimental values are an order of magnitude smaller and in strong disagreement with the 
CHFB+5DCH theory. Thus, for these nuclei at least, the observed band built on the 
states does not correspond to /3-vibration calculated in the CHFB+5DCH theory. 

Recent measurements have shown that many + excited states are present at low excita- 
tion energy in the deformed rare earths Q]. This suggests that the + levels described by 
the CHFB+5DCH may be quite fragmented. For example, the /3-vibrational mode couples 
to such modes as pairing vibrations and/or incoherent 2qp excitations. In the algebraic 
models, efforts have been made to explain the extra + states by introducing many-particle 
many-hole excitations Other regions of deformed nuclei like the actinides and trans- 
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actinides would be worth investigating to check whether they also are missing the coherent 
/3-vibrational structure predicted by the CHFB+5DCH theory. 



3. 2 ~ * Oi transition 

Another observable relevant to the structure of the 2 level is its monopole transition 
strength to the ground state. The strength is conventionally expressed in terms of the 
quantity p 2 defined as [3] 



p 2 (£0;0+ - 0+) = | ( ^ IE g 2 ir?|0t) l 2 , (29) 



•o 



with R = 1.2 A 1 / 3 fm. We calculate the required matrix element as 
z z 
(0^lX>* 2 |°^= / da da 2 g 1 (a ,a 2 )g° %a ,a 2 )ma ,a 2 )\J2r-ma ,a 2 )). (30) 

8=1 J 1=1 

We find rather interesting systematics with respect to the neutron number as shown in Fig. 
[251 The calculations display oscillatory structures with broad maxima located near mid-shell 
closures (N ~ 40, 64, 100, and 150) and sharp minima in the vicinity of shell closures with 
N ~ 20, 28, 50, 82, 126, and 184. Except for light nuclei, all these minima take place for 
mean ground state deformations < (3 > with small values, that is for spherical equilibrium 
shapes. These features are globally consistent with IBA calculations for E0 transitions: the 
p 2 (E0; 2 — > Of) values raise sharply in the shape transition regions and then remain large 
for well deformed nuclei [3]. Other minima in E0 strengths are found below N ~ 126 and 
N ~ 184 neutron shell closures. Those associated with N ~ 162 and weak < (3 > values are 
correlated with the opening of the neutron shell closure in transactinide nuclei, as discussed 
previously in Sec. IIV CI The other two minima take place in open shell nuclei at mean 
ground state deformation in the ranges < (3 > ~ 0.19-0.26 and < (3 > ~ 0.26-0.30 for 
neutron numbers N ~ 116 and N ~ 158, respectively. These features as well as similar ones 
identified in light nuclei with N ~ 20 and N ~ 28 indicate that the existence of minima in 
the E0 strengths over the (N,Z) plane are not exclusively correlated with shell closure and 
spherical ground states. Finally we note that the E0 strength values are at a maximum near 
N = 132 and N = 190 and decrease with N increasing. These features are relevant to the 
Z > 96 isotopic chains for which mean ground state deformations display strong variations 
(see right-hand panel in Fig. [3|). 
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FIG. 25: (Color online) Systematics of the calculated squared monopole transition matrix element 
P 2 (Eq; O2 — > 0~t) as a function of neutron number N. Color code is for proton number. 

We have compared our p 2 (E0) transition strengths to 87 of the 91 nuclei tabulated in Ref. 



771 ] . The result for R figure of merit is given in Table UVl We see that experimental matrix 



elements are on the average very small compared to theory. This suggests that the experi- 
mental levels may have a very different structure than the calculated ones. It may be that 
configurations ignored by the CHFB+5DCH, such as 2 qp excitations, may be important in 



the non-yrast spectrum 



781 ] . This problem with the parameter-free CHFB+5DCH theory 



which globally overestimates the EO strengths by an order of magnitude is shared by other 
models of nuclear structure. For example, enforcing realistic model descriptions of Ml and 
E2 transitions or charge radii, it is not uncommon that calculated EO strengths are up to 
ten times stronger than experimental values 
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80 



8l( | . The actual nature of E0 transitions 



remains an elusive issue pointing to major improvements required in structure models. 
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4- Coexistence between bands 



We have seen above that the conditions for the occurrence of /9-vibration impose the 
medium- and heavy-mass deformed nuclei to lay in specific (Z,N) regions (see Fig. 1241) . 
One is then left with the issue as to what can be learned on collective excited K = band 
properties for nuclei which do not belong to this sample. For this purpose we define an 
important indicator of band structure through the ratio 

R20{BB) - B(E2;2t^0f) (31) 
for all nuclei of present interest with the provision that 2jj~ levels have preponderant K = 
component (i.e. P(K = 0) > 0.75) in their wave functions which unambiguously makes them 
members of excited K = bands. This ratio displays marked structures only if plotted versus 
neutron number. It is shown in Fig. [26] where open and solid symbols are for R^ 2 > 2.3 and 
i?42 < 2.3, respectively. R 2 q{BB') values in the vicinity of R 2 q(BB') =1 are representative of 
the points concentrated at the center point of the triangle shown previously in Fig. [23], and 
these take place near mid-neutron-closed shell numbers N = 100, 150. Most R 2 o(BB') , s take 
on values away from unity. Those with R 20 (BB') > 1 are suggestive of stronger collectivity 
present in excited K = band than in ground state band, and the other way around for 
R 20 (BB') < 1 values. 

It is for N < 60 nuclei that the symbols in Fig. [26] show disparate features as, within a 
narrow range of N values, R 2 o(BB') rapidly flips from R 20 (BB') > 1 to R 2 o(BB') < 1. These 
features have not been analyzed in detail but suggest strong shell effects driving nuclei from 
near spherical to deformed or from prolate to oblate, or vice versa. A typical example is that 
offered by the neutron-deficient Kr isotopes which display shape coexistence features, and 



also undergo a shape transition from prolate to oblate [23j. Nuclides with N > 60 display less 
scattered features in their ratios R 2 q(BB') which now form a seemingly regular oscillatory 
trajectory versus N. Extrema are not well localized but undoubtedly they are reasonably 
close to N = 66, 78, 90, 100, 116, 136, 146, 170 and 196 for nuclei with R A2 > 2.3. Detailed 
information on the location of nuclei with such properties over the (N,Z) plane is not yet 
available. At this moment we have checked that the ratio R 2 q(BB') can serve as good 
indicator for the identification of nuclei known to display shape coexistence (i.e. isotopes 
of the Se, Kr, Sr, Zr, Sm, Hg, Pb, and Po elements) associated with the presence of two 
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FIG. 26: (Color online) The ratio R2o(BB'), Eq. (I31j) . as a function of neutron number. Color code 
is for proton number. Data marked with open and solid symbols are for R42 > 2.3 and -R42 < 2.3, 
respectively. 

(or three) minima in potential energy surfaces. Coexistence between collective bands are 
calculated for Pd, Cd, and Te isotopes which are known experimentally to display such 



features, see e.g. [82j. For these nuclides, coexistence is related to well localized maxima 
in collective masses present at different loci over the (/3,7) plane and not to prominent 
minima in potential energy surfaces. Other theoretical treatments of this kind of coexistence 



have mostly been based on the Interacting Boson Mode 



between shells to explain the intruder states 83|, [8J, l85 |. 



invoking particle-hole excitations 



VII. SUMMARY AND OUTLOOK 



In this work we have evaluated the performance of the CHFB+5DCH theory based on the 
Gogny D1S interaction as a global theory of nuclear structure. Highlights of the successes 
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of the theory are its accurate predictions with respect to charge radii, the classification of 
nuclei as deformed rotors or not, and the ground state band properties of the axial rotors. 

The calculated 2-nucleon separation energies are interesting in that they show shell effects 
and how they are modified for nuclei far from stability. An example is the nuclei near 
N, Z = 40. These are often predicted to be spherical in HFB. This is also the case for the 
Gogny interaction, but we find that there is a change in structure going to the CHFB+5DCH 
extension, and the nuclei become deformed rotors, in agreement with experiment. The 
CHFB+5DCH theory is a suitable framework where the issues of shell erosion and shell 
collapse can be addressed. 

As a spectroscopic theory, the CHFB+5DCH has considerable predictive power for the 
lowest yrast and yrare states, even for nonrotational spectra. This is illustrated by predicted 
energies of the 2^,4+ and l\ excitations, which we compared with compiled empirical data. 
The performance with respect to the levels also showed predictive power, but here sys- 
tematic deficiencies of the theory become apparent. All in all, it is quite remarkable that 
a theory based on a many-body Hamiltonian with only 14 interaction parameters has such 
predictive power over the broad range of the nuclei that can be treated in the methodology. 
An important finding is that deformation alone is not a good predictor of rotational spectra. 
We defined a quantity, /3-softness, that correlates much better. For convenience, we summa- 
rize in Table HVl the figures of merit for the performance of the theory for excitation energies 
and transition properties. We also provide in a retrievable form the specific predictions for 
spectral properties of about 1700 even-even nuclei and analyzed some of the systematics of 
the predicted quantities. 

There are a number of avenues that could be pursued to improve the theory, some of 
which are quite straightforward, at least in principle. The treatment of the inertial masses 
could be improved by using the Thouless-Valatin prescription which is better justified than 
the cranking approximation we have used up to now. This would surely improve the energies 
of the 2 excitations, which is one of the problems of the current theory. This would require 
calculating the QRPA response function at every grid point. This would of course add to the 
computational burden, but the ingredients to perform the calculation are available for the 
most part in the intermediate calculations already performed. The rotational inertias could 
also be improved by using a finite amplitude rotational field uJk adjusted to self-consistency 
for the angular momentum value being calculated 86|. The ingredients of this self-consistent 
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TABLE IV: Summary of the performance statistics of the CHFB+5DCH for excitation energies 
and transition properties. See text for definitions of R and a. The column "Number" gives the 
number of nuclei in the comparison data set. 



Observable 


Number 


R a 


E(2t) 


513 


0.11 0.35 


B(E2;2f^0t) 


311 


0.20 0.42 


R42 


480 


0.03 0.14 


R&2 


427 


0.08 0.21 




352 


0.19 0.30 


E(0t) 


317 


0.31 0.36 




87 


2.1 1.9 



treatment are already available and have been implemented in some calculations [171 ] . This 
also would considerably add to the computational burden. 

We have taken the correlation energy as an indicator of the validity of the Gaussian 
Overlap Approximation, and excluded nuclei whose calculated correlation energies are un- 
physical. A better treatment of the correlation energy would include the ZPE term coming 
from the curvature of the potential energy surface. This term is small in nuclei with broader 
collective potential energy surfaces, but in noncollective nuclei the minima can be narrow, 
and this potential term in the Hamiltonian mapping might have a significant effect. 

Another deficiency of the CHFB+5DCH that became evident in the discussion of the 
level properties is the need for 2 qp components in the wave functions. This can be carried 
out in the GCM if the Hamiltonian operator in the collective space is calculated with full 
treatment of the nonlocality js7]. However, there is no clear road to us for how to include 
these components and keep the Gaussian Overlap Approximation as yet. 

While it is remarkable that the Hamiltonian based on the Gogny D1S interaction has 
so much predictive power after 30 years, one can still ask how the calculated observables 
depend on the intrinsic properties of the Hamiltonian and whether an improvement can be 
made at that level. We note that there has already been some work to improve the Gogny 
interaction for calculating nuclear masses, while keeping unaltered the performance of the 
CHFB+5DCH and RPA theories achieved previously with D1S Q|. 
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Finally, as a separate publication, we will present in more detail the spectroscopic prop- 
erties to higher angular momenta for the deformed rare-earth nuclei. This will also include 
discussion on the systematica of interband transitions and odd-even angular momentum 
staggering of the 7 bands 
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